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A study of the fine scale motions of 
incompressible time-developing mixing layers 


By J. Soria 1 , M. S. Chong 2 , R. Sondergaard 3 , 
A* E. Perry 4 AND B. J. Cantwell 3 


This work is an extension of a project conducted at the previous CTR summer 
program and was reported by Chen et al. (1990). In that program, the geometry 
and topology of the dissipating motions in a variety of shear flows was examined. All 
data was produced by direct numerical simulations (DNS). The partial derivatives 
of the velocity field were determined at every grid point in the flow and various 
invariants and related quantities were computed from the velocity gradient tensor. 
Motions characterized by high rates of kinetic energy dissipation and high enstrophy 
were of particular interest. Scatter diagrams of the invariants were mapped out and 
interesting and unexpected patterns were seen. Each type of shear layer produced 
its own characteristic scatter plot. 

In the present project, attention is focused on the incompressible plane mixing 
layer, and the scatter diagrams are replaced with more useful joint probability 
density contours. Comparison of the topology of the dissipating motions of flows 
at different Reynolds numbers are made. Also, plane mixing layers at the same 
Reynolds number but with different initial conditions are compared. 


1. Method of approach 

The velocity gradient tensor may be broken up into a symmetric and an antisym- 
metric part Aij = dutjdxj — Sij + Wij where Sij = (dui/dxj + duj/dxi )/ 2 and 
Wij = ( dui/dxj - dujldx t )/2 are the rate-of-strain and rate-of- rotation tensors, 
respectively. The eigenvalues of Aij satisfy the characteristic equation 

A 3 + PA 2 + Q\ + R = 0 (1) 

where the matrix invariants are: 

P = — (An "t" A 22 “h A 33 ) = trace [A] — Sa (2) 
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It can be shown that, in the P — Q — R space of matrix invariants, the surface 
which divides characteristic equations with three real solutions for the eigenvalues 
from characteristic equations with one real and two complex solutions is 


27 R 2 + (4P 3 - 1 SPQ)R -I- (4<? 3 - P 2 Q 2 ) = 0. (5) 


A detailed discussion of the properties of this surface is given in Chong, Perry & 
Cantwell (1990) along with a guide to the various possible elementary flow patterns 
which can occur in different domains. 

Much of the discussion in this report concerns the symmetric part of the velocity 
gradient tensor, the second invariant of which is proportional to the negative of the 
kinetic energy dissipation. The invariants of the rate-of-strain tensor, P5, Qs and 
i?s, are generated by setting the components of Wij to zero in the above relations. 
The flows considered are, with one exception, incompressible hence P = Ps = 0. 
Thus the local geometry of the flow is completely described by the second and third 
invariants (Q, R) and {Qs, Rs )• The second invariant of the rate-of*rotation tensor, 
Qw > is non-zero and is proportional to the enstrophy. The first and third invariants 
of the rate-of-rotation tensor are identically zero. 

The method for classifying the flow structure was first developed at the 1990 
CTR summer program by Chen et al. (1990) and is described as follows: 

(i) Evaluate the nine partial derivatives of the velocity gradient tensor at every point 
in the computed field. 

(ii) Evaluate Q, ij, Qs, Rs and Qw at every point. 

(iii) Create scatter plots of the results in the space of invariants, Q versus R, Qs 
versus Rs , and —Qs versus Qw • 

Figure 1 illustrates the various flow topologies which can occur in the plane P = 0. 
The intersection of this plane with the surface (5) is given by 
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FIGURE 1. Three-dimensional topologies in the Q — R (P — 0) plane. 

which divides real solutions from complex solutions as indicated. 

For the case P — 0, the second invariant is 

Q = \[W ij W ii -S tj S ij ] (7) 

where the indices have been switched to indicate explicitly that Q is formed from 
the difference of two terms, each of which is a positive sum of squares. The local 
topology has complex or real eigenvalues depending on whether the (Q, R) pair 
evaluated at a given point in the flow lies above or below (6). 

The mechanical dissipation of kinetic energy due to viscous friction is 

<f> = 2 vSijSij = —4 vQs- 


( 8 ) 
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(b) 



(c) 


(d) 


Figure 2. Scatter plots of (a) Q vs. R, (b) Q s vs. R 3 , (c) —Q s vs. Q w for hipairex 
at tU/So = 29.8, and (d) Q 9 vs. Q w for a compressible mixing layer computed by 
Chen (1991) at tU/6 0 = 72.0. 


Large negative values of Qs correspond to large rates of dissipation of kinetic energy. 
Large negative values of Q indicate regions where the strain is both large and 
strongly dominant over the enstrophy. Large positive values of Q indicate the 
reverse. 


2. Results 

We will consider in this paper the incompressible (P = 0) plane mixing layers 
computed by Rogers and Moser at NASA Ames. Three direct numerical simulations 
(DNS) are considered, namely, hipairex , mega , and tbi The cases hipairex and mega 
were initiated from laminar error function profiles, and tbl was initiated with two 
turbulent boundary layer realizations with equal and opposite free stream velocities 
placed on opposite sides of a dividing plate which was dissolved at time t = 0. The 
initial turbulent boundary layers were DNS computations of Spalart (1988). 
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FIGURE 3. Planar x-z average Q 9 vs. cross-stream direction y for hipairex at 
tU/6 0 = 22.3,25.3,29.8. 



FIGURE 4. Weighted Probability Density Functions P( — Q s ) and (— Q S )P(— Q s ) 
vs. —Qs for hipairex at at W/Sq = 29.8. 
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FIGURE 5. Vorticity thickness Reynolds number Rs vs. non-dimensional time for 
the incompressible mixing layers. 

Table I 



hipairex 

mega 

tbl 

A U = U 2 - Ui 

2 

2 

2 

Initial vorticity thickness, 8q 

1 

1 

1.4 

Viscosity, v 

1/250 

1/250 

1/500 

Initial Re = AU8o/2u 

250 

250 

700 


All cases were computed as time developing layers and Table I shows the prop- 
erties of the layers. Details of the hipairex results have been reported by Moser &: 
Rogers (1990) and Rogers &: Moser (1992). Unless otherwise stated, all results are 
normalized by half the velocity difference across the layer, U , and the initial vortic- 
ity thickness 8q. Figures 2(a), 2(b) and 2(c) show scatter diagrams taken from Chen 
et al. (1990). These diagrams are made up of the entire data set for a given time. 
Figure 2(a) shows that most of the high gradient motions belong to the topology of 
stable focus stretching. Figure 2(c) is most informative. Data which falls on the line 
of 45° through the origin represents high dissipation accompanied by high enstro- 
phy. It can be shown that such points come from vortex sheets where most of the 
rate-of-strain is dominated by the velocity gradient within the sheet. Data which 
lies along the horizontal axis represents high enstrophy with little dissipation as 
would occur in solid body rotation in vortex tubes. As a matter of interest, Figure 
2(d) shows a plot from a compressible mixing layer computation by Chen (1990), 
and, according to the figure, the data could be described as primarily sheet-like. 
The reason for this is a mystery at this stage. 




108 J. Soria , M, S . Chong, R. Sondergaard, A. E. Perry, & B. J. Cantwell 


Qs Qs 



Figure 7. Q 3 vs R a for hipairex. tU/6 = (a) 19.3, (b) 22.3, (c) 25.3, (d) 29.8. 


Classical arguments, based on the idea that dissipation of turbulent kinetic energy 
scales with production, lead to the following estimates: 


« - 


(9) 


where the 5 t -j are fluctuating non-normalized strain rates. Results from experiment 
show that for fully developed shear layers 


-uV/(AC7) 2 ~ .012 = 


4 U 2 ' 


( 10 ) 


From (9) and (10) it can be shown that 


e = 2 vS'ijS'ij = meu 3 /s 


( 11 ) 
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FIGURE 8. — Q, vs Q w for hipairex. tU/6 = (a) 19.3, (b) 22.3, (c) 25.3, (d) 29.8. 

For time tU/6 0 = 29.8 in hipairex where the vorticity thickness has increased by a 
factor of 6.5 over the initial thickness, the Reynolds number Rg based on the current 
vorticity thickness, 6, and the velocity difference across the layer is 3000. Hence 

(M 2 

U 2 ~ 4 6 \6 ) 


= 1.704. 


( 12 ) 



110 J. Soria , M. S . Chong , R. Sondergaard, A . E. Perry , & B. J. Cantwell 



Figure 9. Q vs R for me^a. tU/6 = (a) 21.0, (b) 25.0, (c) 35.0, (d) 49.0. 

One would expect the average value of this quantity at the midplane of the mixing 
layer to be of this order. 

Mean profiles of —Qs at various times are shown in Figure 3 and are half the value 
indicated by (12) which is indicative of the production of kinetic energy. This ratio of 
about 2 for production to dissipation is in agreement with the fully developed value 
obtained from experiments by Bradshaw & Ferriss (1967). An order of magnitude 
analysis similar to (12) reported by Chen et al. (1990) giving the value of 18.2 was 
in error due to incorrect normalization of the variables. 

Figure 4 shows the weighted probability density function of — Qs over the entire 
volume of the mixing layer, and most of the contribution comes from — Qs between 
0 and 3. Although the far flung values of — Qs on the scatter diagram tend to 
follow interesting patterns, they contribute only of order 10% to the total energy 
dissipation. For this reason, it was felt that scatter diagrams should be replaced 
by joint probability density diagrams with contours corresponding to the logarithm 
of the probability density function so that possible ridges could be seen in regions 
which are highly darkened in the scatter plots. 
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FIGURE 10. Q » vs R„ for mega. tU/6 = (a) 21.0, (b) 25.0, (c) 35.0, (d) 49.0. 

From (11) and (12) it can be seen that Qs normalized by the current vorticity 
thickness scales with Re, and, therefore, it seems likely that Rs should scale with 
R 3 / 2 . This would imply that the data should follow a curve 

|b.1«(I0.I) ! (13) 

This relationship is what one might expect purely on dimensional grounds, but 
there is no rigorous proof. It is interesting to note that such a curve on the Qs 
versus Rs plot represents a rate of strain geometry where the principal rates of 
strain a, (3 and 7 are in a constant ratio to one another. For the data set hipairex, 
points of high dissipation follow closely the curve corresponding to the ratio of 
a: /?: 7 = 3:1: —4, which was observed by Ashurst et al. (1987) in studies of 
forced isotropic turbulence. In addition, as noted by Sondergaard et al. (1991), the 
vorticity vector tends to align itself with the second principle rate of strain (3. It 
should be noted that while other data sets analyzed by Sondergaard et al. (1991) 
show the same vorticity alignment, the 3: 1: —4 ratio of rates of strain is not always 
observed. 

The result depicted in 2(b) is that motions characterized by very high rates of 
dissipation (large negative Qs) clearly show a preference for the right half plane of 
Figure 2(b) corresponding to a local topology of the rate of strain tensor which is 



112 J. Soria , M. S. Chong, R. Sondergaard, A. E . Perry, & B. J. Cantwell 



FIGURE 11. — Q 3 vs Q w for mega . = (a) 21.0, (b) 25.0, (c) 35.0, (d) 49.0. 

of the type saddle-saddle-unstable node (cf. Figure 1). From Figure 2(a), it can be 
seen that the velocity gradient tensor admits all possible incompressible topologies 
although there is, nevertheless, a great deal of structure in Figure 2(a). Not only 
is the basic scaling (12) observed, but it appears that, with a modest amount of 
scatter, the fine scale motions follow a relation of the form 

R 9 ^K(-Q 9 )$. (13) 

The positive quantity K is expected to be a function of the Reynolds number with 
an upper limit of K = 2>/3/9 corresponding to locally axisymmetric flow (cf. Figure 
i). 

3. Comparison with high Reynolds number flows 

Figure 5 shows a plot of Reynolds number based on current vorticity thickness 
for the three cases mentioned earlier. Figures 6, 7, and 8 show the invariant plots 
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Figure 12. Q vs R for tbl. tU/6 = (a) 0.0, (b) 26.25, (c) 47.5, (d) 76.25. 


for hipairex in contour plot form for the joint probability density distributions. 
Figures 9, 10, and 11 show the results for mega. These results are similar to the 
scatter diagrams given in Figure 2 but are corrected for nonuniform grid spacing. 
More structural features are apparent, and an interesting feature in Figures 6 to 8 
is that velocity gradients tend to increase with time and at the latest time show 
a decrease. In fully developed turbulent plane mixing layers, if dissipation scales 
with production, then according to the Kolmogorov scaling, the velocity gradients 
should decrease with time. According to this reasoning, hipairex is under-developed 
for most if not all of the times shown. It is unclear whether, at the latest time, the 
gradients are beginning to decrease because the flow is reaching a fully developed 
state or because of constraining by the grid. Figures 9, 10, and 11 show the results 
for mega , and there appears less pronounced sheet-like structures but more tube-like 
patterns for the higher Reynolds numbers. 

Figures 12, 13, and 14 show similar results for tbl. which started out as two 
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Figure 13. Qs vs J? 4 for tbl. tU/6 = (a) 0.0, (b) 26.25, (c) 47.5, (d) 76.25. 

turbulent boundary layers and then developed to a much higher Reynolds number 
than hipairex. Figure 12 is most interesting. It shows that all data points for 
the turbulent boundary layer cluster near the origin of the Q versus R plot and 
suddenly explode to much higher gradients in the plane mixing layer. These pictures 
graphically illustrate how much greater velocity gradients become when the wall 
constraint of a turbulent boundary layer is removed. It should be noted that near the 
wall, the Q and R of a turbulent boundary layer are small even when the gradients 
aren’t. A better measure of the relative magnitudes of the velocity gradients can 
be inferred from figure 14. Again, the gradients tend to grow and then diminish 
at late times. The Qs versus Rs plot shows that the strain rates tend to follow a 
different curve, closer to the real-imaginary dividing surface (6). Hence, this aspect 
of the fine scale motion appears to be Reynolds number dependent. The plots in 
Figure 14 show that the turbulent boundary layer structures at t* = 0 are sheet-like, 
but, in contrast to hipairex , there are no preferred structures revealed by the — Qs 
versus Qw plot for later times (cf. Figure 8). As with hipairex the Q versus R plot 
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FIGURE 15. Vortex lines for tbl at tU/8 = 0.0. 
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FIGURE 16. Vortex lines for tbl at tU/6 = 76.3. 



FIGURE 17. Streamlines for tbl at tU/6 = 76.3. 

shows that most of the gradients belong to the topology of stable focus stretching. 

The highly organized patterns seen in Chen et al. (1990) for hipairex are replaced 
by most complex structures in tbl. Vortex lines for tbl are shown in Figures 15 and 
16. Shown in Figure 15 is the initial turbulent boundary layer, and the attached 
eddies which lean approximately 45° to the mean flow direction are apparent. In 
Figure 16 are shown vortex lines of the plane mixing layer after some development. 
Although no clear spanwise rollups are apparent from this vorticity plot, Figure 17 
shows instantaneous streamline patterns which indicate possible large scale spanwise 
roll-ups. 

4* Comparison of two initial conditions at the same Reynolds number 

From Figure 5, it can be seen that there is an overlap of Reynolds numbers for 
mega and tbl. In fact, they both share a Reynolds number of 5000 as indicated in the 
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Figure 20. Rescaled ~Q S vs. Q w plots for (a) mega at tU/8 = 49.0, (b) tbl at 
tU/8 = 61.0. 

figure. Figures 18(a) and (b) show Q versus R plots of tbl and mega , each scaled with 
the current vorticity thickness and appropriate velocity U. The joint probability 
density contours have been rescaled to account for the different population density 
of points. Figures 19 (a) and (b) show the corresponding -Qs versus Qw plots for 
comparison. Although the shape of the plots are roughly the same, there appears 
to be a major difference in the scaling, indicating that the velocity gradients in tbl 
are considerably lower than in mega for the same Reynolds number. The reasons 
for this difference need to be pursued in future work. 

5, Conclusions 

In all flow cases considered here, motions with the highest dissipation of ki- 
netic energy per unit volume were of the topological classification stable focus with 
stretching as found from Q — R plots. 

In the case designated as hipairex, the flow was initiated from a laminar layer with 
an error function profile and the maximum Reynolds number Rs to which the flow 
evolved was 3000. Here, the highly dissipative motions were usually accompanied 
by a high enstrophy indicating a vortex sheet-like structure. From Qs versus Rs 
plots, the rate of stain tensor for dissipating points had a topology of unstable node 
saddle-saddle with the rate of strains being of a given ratio and with the vorticity 
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vector tending to align with the intermediate strain /?. 

In the case of tbl , which was initiated from two turbulent boundary layers placed 
back to back, the highest local Reynolds number considered was Re — 9000. Here, 
the —Qs versus Qw plots indicated no preferred structure for the highly dissipating 
motions although the Q — R plots indicated a strong preference for stable focus 
stretching. Also, the Qs versus Rs plots showed that the highly dissipating motions 
tend toward a: (3: y = 1 : 1 : —2. No vorticity alignment checks were made, but it is 
expected that the vorticity vectors will tend to align with the (3 axis (cf. Sondergaard 
et a/., 1991). 

Comparison of two flows at the same local Reynolds number but with two entirely 
different initial conditions was made using flow cases designated mega and tbl . Plots 
of Q versus R and — Qs versus Qw when nondimensionalized appropriately show 
essentially the same topological structure and scaling from Re = 5000 even though 
mega was initiated from a laminar error function profile layer and tbl from turbulent 
boundary layers. Although the shape of the plots are roughly the same, there 
appears to be a major difference in the scaling, indicating that the velocity gradients 
in tbl are considerably lower than in mega for the same Reynolds number. The 
reasons for this difference need to be pursued in future work. 
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